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ABSTRACT 

Results  of  calculations  are  given  for  lateral  drifting  forces  acting 
on  a cylinder  and  a Series  60  hull  derived  by  a new  procedure  involving 
application  of  the  quadratic  frequency  response  function  (QFRF)  and  the 
close-fit  method  for  flows  inducted  by  hull  sections  in  the  near  field.  The 
near  field  inductions  are  required  by  the  QFRF  in  order  to  obtain  the  non- 
linear (second  order)  interaction  effects  in  the  presence  of  dual  waves. 
This  method  yields  the  mean  drifting  force  consisting  of  four  components  of 
which  the  relative-wave-elevation  term  is  dominant,  whereas  the  Bernoulli 
quadratic  term  is  secondary.  A mathematical  discontinuity  of  the  drifting 
force  in  the  neighborhood  of  a high  frequency  is  investigated  by  applying  a 
modified  Green's  function.  Mean  drifting  forces  for  two  forms  have  been 
calculated  and  the  results  compared  with  available  analysis  and  model  data. 
A fairly  complete  quadratic  frequency  response  function  for  lateral  drift- 
ing forces  in  the  bi-frequency  domain  was  computed  for  the  Series  60  Ship 
Cg  = 0.6  at  sixty-degree  heading.  The  computed  results  are  presented  in  a 
three-dimensional  view  accompanied  by  a detailed  discussion  of  the  charac- 
teristic behavior. 
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NOMENCLATURE 


Va2 


wave  amplitude 


amplitude  of  wave  component  of  frequency  uu , w u)  , w 

I 2 m n 


A_  wave  amplitude  ratio 

B beam  of  a section  in  the  load  waterplane 

c(x)  a transverse  section  contour 

Cp  non-dimensional  lateral  force  coeff ic iert  (Eq.26) 

Cu  non-dimensional  lateral  moment  coefficient  (Eq.27) 

O^Ct)  drifting  force 

e wave  elevation  (Eq.20) 

f(x)  sectional  lateral  force 

f mean  sectional  lateral  force 

F lateral  force 

F mean  lateral  force 

9j(t)  linear  impulse  response 

g2(tr t2)  second  degree  kernel  analogous  to  an  impulse  response 

G pulsating  source  function 

Gj(ai)  linear  frequency  response  function 

^2^m,ia>n^  quadratic  frequency  response  function  defined  by  Eq.(3^) 

h wave  height  or  wave  elevation 

hi  t-(a)  i^  ) functions  defined  by  Eq.(32) 

I , k m n 


1^  integral  expression  defined  by  Eq.(3l) 

L hull  length 

mQ(x)  mass  of  a section 

M lateral  moment 

M mean  lateral  moment 
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pressure 


source  intensity 
relative  wave  elevation 


displacement  of  a response  motion  of  mode  m 
segment 


T draft,  period 

x,y,z  moving  body  coordinates 

X,Y,Z  space  fixed  coordinates 

A displacement 


wave  number 


wave  number  of  frequency 
water  dens i ty 

resultant  velocity  potential 
diffraction  potential  of  mode  m 


radiation  potential  of  mode  m 


"•W 


circular  frequency 


Q.  circular  frequency 

a non-dimensional  frequency 


non-dimensional  time 
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INTRODUCTION 

The  ultimate  goal  of  the  present  study  is  to  establish  an  improved 
prediction  technique  of  the  fluctuating  lateral  drifting  force  and  moment 
(or  second-order  yaw  moment)  and  then  to  determine  their  statistical  charac- 
teristics. Thus,  a useful  tool  would  be  available  for  prediction  of  the 
behavior  of  naval  vessels  in  random  seas. 

The  broad  problem  is  divided  into  two  Phases:  Phase  I provides  an 
analysis  and  limited  correlation  with  model  measurements  of  mean  drifting 
force;  Phase  II  was  proposed  to  be  directed  to  a complete  experimental 
determinat ion  of  the  QFRF  and  subsequent  correlation  with  the  theoretical 
predictions  evolved  in  Phase  I.  This  report  deals  only  with  Phase  I. 

A recent  investigation*  of  the  quadratic  frequency  response  for  added 

resistance  has  shown  satisfactory  agreement  with  the  results  of  an  experi- 

2 

mental  investigation.  Reference  2 shows  that  the  quadratic  response  function 
can  be  utilized  to  predict  reliably  the  very  low-frequency,  fluctuating  added 

resistance  in  irregular  waves.  The  analytical  technique  developed  in  these 

I 2 3 

studies  ’ has  been  applied  by  Kim  and  Breslin  for  prediction  of  slow  drift 

oscillation  of  a moored  ship  in  random  head  seas.  This  method  will  be  employed 

in  the  present  investigation. 

U 

Simplified  analytical  methods  have  been  proposed  by  Hsu  and  Blenkarn, 
Newman,'’  and  Faltinsen  and  Loken,^  for  the  prediction  of  slow  drift  oscilla- 
tion. These  procedures  depend  on  the  mean  drifting  force,  which  represents 
only  a special  output  from  ORFR  procedures. 

Maruo^  provided  a basis  for  predicting  the  mean  drifting  force  due  to 

g 

a monochromatic  wave  on  a cylinder  floating  in  a beam  wave  train.  Ogawa 

applied  Maruo's  method  to  a fixed  ship  by  using  a strip  method  based  on  the 

concept  of  the  snake-type  wave  generator.  He  showed  that  his  calculation 

q 

agreed  well  with  his  experiments.  Spcns  and  Lalangas  reported  experimental 
results  on  the  mean  lateral  drifting  force  and  moment  on  a ship  model  in 


Superior  numbers  in  text  matter  refer  to  similarly  numbered  references  listed 
at  the  end  of  this  report. 
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oblique  waves.  Newman  calculated  the  mean  lateral  drifting  force  based 

on  the  slender  body  theory.  However,  his  predictions  did  not  agree  well 

9 1 1 

with  the  experiments.  Kim,  et  al,  calculated  the  mean  lateral  drifting 

force  and  moment  on  a fixed  as  well  as  a free  ship  in  oblique  waves,  using 

8 12 

a strip  method  and  showed  acceptable  comparisons  with  experiments,  ’ 

13 

Salvesen  developed  a method  of  calculating  the  mean  added  resistance  and 

mean  lateral  drifting  force  by  using  a strip  method.  Compared  with  the  re- 

9 

suits  of  Spens  and  Lalangas,  his  results  showed  an  appreciable  discrepancy 
He  attributed  this  discrepancy  to  difficulties  encountered  in  the  experiment 
and  to  "weak  scatter"  assumptions  made  in  Newman's  theory.'^ 

1 4 

Pinkster  and  Oortmerssen  applied  a three-dimensional  source  dis- 
tribution method  to  the  evaluation  of  the  mean  drift  forces,  which  were 
obtained  by  integrating  the  pressure  distribution  on  the  wetted  hull  surface. 

The  present  calculations  are  carried  out  in  two  steps;  (A)  calculation 
of  the  lateral  second-order  force  acting  on  a ship  section,  free  or  fixed  in 
regular  beam  waves,  and  (B)  stripwise  calculation  of  the  quadratic  frequency 
response  function  for  the  ship  in  oblique  dual  waves.  The  analysis  is  based 
on  the  hypothesis  that  the  contribution  of  the  second-order  potential  to  the 
second-order  force  is  negligible  compared  to  those  contributed  by  the  first- 
order  potential.  Rigorous  analyses  of  the  second-order  force  and  moment  on 

a cylinder  free  in  waves  have  been  given  by  Potash^  and  Soding,^  both  of 

£ 

which  do  not  neglect  the  second-order  potential.  Faltinsen  and  Loken  also 
estimated  the  second-order  potential  in  their  evaluations  of  the  slow  drift- 
ing force  and  moment. 

The  present  analysis  integrates  the  pressure  distribution  on  the 
section  surface  (as  in  Reference  14),  whereas  Maruo^  integrates  pressures 
over  control  surface  in  the  far  field.  The  former  method  may  therefore  be 
designated  as  a "near-field  method, " wh i le  the  latter,  a "far-field  method." 

The  far-field  method^  yields  a simple  formula  for  the  mean  force  whereas 
the  near-field  method  presents  a sum  of  several  components.  The  latter  seems 
to  be  complicated  at  first  and  difficult  to  use;  however,  it  clearly  reveals 
the  effects  of  the  behavior  of  the  response  motions  such  as  heave  and  roll 
resonances  and  body  parameters.  Effects  of  motions  and  body  parameters  have 
remained  unknown  and  sought  after  for  many  years  by  researchers.  Moreover, 
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the  far-field  method  cannot  be  used  for  the  analysis  of  the  quadratic  fre- 
quency response  for  the  drifting  force  and  moment. 

A two-dimens ional  procedure,  mentis  ad  above,  has  been  employed  as  in 
Reference  11,  for  stripwise  calculation  of  the  force  and  moment  of  a ship  in 
oblique  dual  waves. 

A special  output  of  the  quadratic  frequency  procedure  is,  in  fact,  the 
mean  drift  force  or  moment  as  demonstrated  later.  Numerical  calculations  of 

the  mean  drifting  force  and  moment  have  been  carried  out  for  a cylindrical 

17  12 

model,  tested  by  Koterayama,  and  for  a ship  Series  60,  C =0.6  by  Lalangas. 

B 

The  present  analysis  has  been  compared  with  the  calculation  by  Maruo's  method 
and  model  test  data.  Excellent  agreement  between  the  two  methods  has  been 
confirmed  for  the  bodies  in  beam  seas.  In  oblique  seas,  however,  the  two 
analytical  methods  appear  to  be  in  good  agreement. 

The  mean  drifting  force  or  moment  consists  of  four  components,  i.e., 
gy roscop i c-coupl i ng,[ force  rate]  x [displacement],  the  Bernoulli's  quadratic, 
and  the  relative  wave  terms. 

A highlight  of  the  analysis  of  the  components  shows  that  the  relative 

wave-induced  drifting  force  component  is  the  most  dominant,  seconded  on  ly  by 

Bernoulli's  quadratic  pressure- i nduced  term.  In  addition,  the  influence  of 

the  peak  hull  displacements  on  the  component  forces  have  been  clearly  observed. 

The  presence  of  a mathematical  discontinuity  at  a high  frequency  has  been  in- 

1 8 

vestigated  by  employing  Og i 1 v i e-Sh i n ' s modified  Green's  function. 

The  quadratic  frequency  response  functions  for  drifting  force  and  moment 
of  ship  Series  60,  Cg=0.6  have  been  calculated  at  a sixty-degree  heading  and 
the  results  are  illustrated  in  a three-dimensional  view  (Figure  6). 

The  successful  application  of  the  quadratic  frequency  procedure  which 
yields  the  mean  drifting  forces  (and  moments)  as  a specialized  output  lends 
support  to  the  validity  of  the  technique.  Further  confirmation  of  the  utility 
of  this  procedure  for  elucidation  of  the  bi-frequency  domain  characteristics 
of  the  entire  QFRF  (as  described  in  Phase  II  in  Davidson  Laboratory  Proposal 
P-17I3R)  requires  model  measurements  currently  unavailable. 
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THE  INPUT-OUTPUT  MODEL 


A few  definitions  and  relations  pertaining  to  the  general  input- 
output  rode  1 are  needed  in  the  present  work.  The^e  will  be  given  here  for 
convenience,  although  they  may  be  found  in  Reference  2 in  greater  detail. 

The  general  input-output  model  is  a time  domain  model  and  is  written 


as : 


Fw(t)  = J 9;  ( t, , )'H(t-t1  )dt. 

+ II  g2(ti*ts)Tl(t-t1)Tl(t-ta)dt1dt; 


(1) 


(Omission  of  limits  on  the  integrals  signifies  limits  of  £»  ,) 

In  Eq.  (1):  F (t)  is  the  force  or  moment  exerted  by  the  waves 
w 

Tl(t)  is  a suitably  chosen  wave  elevation  (a  zero  mean  input) 

9,  (t)  is  a 1 inear  impulse  response 

9_  (t,  is  a second  degree  kernel  analogous  to  an  impulse 

response 

With  respect  to  application  to  drifting  force  and  moment,  the  ship  is  assumed 
to  be  station-keeping  in  a wave  system  (the  input)  which  is  defined  as  the 
wave  elevation  at  a point  stationary  with  respect  to  the  mean  position  of  the 
ship.  The  second  degree  kernel  is  assumed  to  be  symmetric  in  its  arguments, 
and  both  kernels  are  assumed  smooth  and  absolutely  integrable. 

The  first  term  in  Eq.(l)  is  the  ordinary  linear  convolution  of  a linear 
impulse  response  function  g^tj)  with  T)(t).  The  second  term  is  a double  con- 
volution, and  the  second  degree  kernel  92{tj,ta)  might  be  called  a quadratic 
impulse  response.  According  to  basic  assumptions,  the  two  impulse  responses 
are  related  to  the  linear  and  quadratic  frequency  response  functions  by 
Fourier  transforms.  These  two  sets  of  transform  pairs  may  be  defined  as 
fol lows ; 


t 

i'.- 


r 
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91  (T)  = I C+IUn  Gi  (u,)d|X 

Gi  (“)  = J e 11,1  9j  (t) dT  (2) 


9?(v’t?)  = JJ'Expi>itvl+iv:> 


Gs^U5i,U?)  = JjExp[- iui  t,  - iu).,T_  ]g?  (Tj 


]G.  (u^  , a'_ ) du;  dir„ 
,T^)di  dt^ 

c * r 


(3) 


The  linear  frequency  response  function,  Gj(a!),  defined  by  Eq.  (2) 
is  absolutely  conventional.  The  quadratic  frequency  response  function, 
G2(u.'j,U).J  is  defined  in  a b i - f requency  plane.  Because  the  kernel  g 
is  assumed  to  be  symmetrical  in  its  arguments: 

9S  (Tj  >t2)  = 92(VS)  (*0 

and  thus: 

G_  (“ij  >0  = G_  (o>_  ,u^  ) (5) 

= G (-U1  ,-«!„) 

e • c <■  J c 

= (6) 

(The  asterisk  denotes  the  complex  conjugate.) 

Figure  la  summarizes  the  results  of  applying  Eq.s  (5)  and  (6)  in 
the  bi-frequency  plane  for  the  eight  possible  coordinate  positions  of  two 
frequencies  whose  absolute  values  are  a and  b.  Equation  (5)  results  in  a 
line  of  synmetry  along  the  line  -u^  . Equation  (6)  results  in  a line  of 
symmetry  of  the  real  part  of  G^u^,^)  defined  by  u)?"-iCj  (and  it  may  be 
noted  that  along  this  line  the  imaginary  part  of  the  function  is  zero). 
These  two  lines  and  the  , u:_  axes  divide  the  bi -frequency  plane  into 
octants,  of  which  the  two  on  either  side  of  the  positive  u.^  axis  may  be 
arbitrarily  chosen  for  reference.  The  equalities  shown  in  Figure  la  in 
the  remaining  six  octants  are  arrived  at  by  applying  Eqs.  (5)  and/or  (6). 
The  assumptions  of  synmetry  of  the  second  degree  kernel  result,  with 
Eq.  (3),  in  a complete  definition  of  G (m  , u>  ) if  the  functions  are  defined 
in  any  pair  of  octants  including  a semi-axis  of  either  frequency.  Thus 


k?' 


i 


S 
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without  loss  in  generality  interpretation  of  the  quadratic  frequency 
response  needs  only  to  involve  the  octants  on  either  side  of  the  positive 
<jUj  axis.  In  these  octants  U):  is  positive  and  | u.^  I ? j | . 

The  interpretation  of  the  quadratic  frequency  response  function  is 
less  direct  than  for  the  linear  case,  but  can  be  approached  in  a similar 
manner.  If  in  the  linear  case  the  system  is  considered  to  be  excited  by 

T| ( t ) = a cosiut 

the  output  may  then  be  written: 

RejaGj  (u>)  Exp(  i tut ) | 

and  Gj  (u>)  is  interpreted  in  terms  of  normalized  amplitude  and  phase  of 
response. 

To  interpret  the  quadratic  frequency  response,  dual  harmonic  exci- 
tation is  necessary.  Accordingly,  it  may  be  assumed  that: 

T|(t)  = a,  cosui  t + a_cosu:_t  (7) 

• * c c 

In  accordance  with  the  previous  discussion  of  symmetry,  both  frequencies 

) are  considered  positive  and  | j s | U.'?  | . The  basic  model,  Eq.  (1) 
is  good  for  any  zero-mean  excitation.  Accordingly,  Eq . (7)  may  be  substi- 
tuted directly  in  Eq . (1).  After  some  algebra  the  final  result  for  the 
response  to  dual  harmonic  excitation  may  be  written  as  follows: 

F?  3 Rclai  Gi  ) E*P(  t)+as  Gj  (tr?)  Exp(  iu^  t)  j 

+ ^ Re{a‘ G_  (a^  ) Exp(  iSuij  t)  | 

+ ^ Re{a^Gs  (a-„,u;?)  Exp(  i2a;?t)  ] 

+ Rejaj  asGg  (^  , ) Exp[  i +u>8)  t] j 

+ Reja,  a?Gs  (aij  ,-a's)  Exp[  i (a^ -u,_)  t]|  (8) 

This  result  shows  that  the  response  of  the  quadratic  system,  Eq . (1), 
to  dual  exc;tation  contains,  in  general,  a shift  in  the  mean  and  components 
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of  six  different  frequencies  [u^  , a'?,  2u)^ , 2c u?,  (u!j+u^),  and  (u^ -U)  )]. 

The  first  two  terms  of  the  result  are  the  superposition  of  the  linear 
responses  at  the  excitation  frequencies.  The  third  and  fourth  terms  of 
Eq.  (8)  represent  a shift  in  the  mean.  These  terms  allow  the  identification 
of  the  mean  drifting  force  or  moment  operator  as  the  value  of  Gs^.Wa)  along 

the  line  u>s=-u;  (or  G„  (u^ , — uUj ) ) . The  fifth  and  sixth  terms  are  the  second 
harmonic  components  (211^,2(1^).  Similarly,  these  terms  allow  the  identifi- 
cation of  second  harmonic  response  with  the  values  of  G^u^jUJ  ) along  the 
•ine  (or  G2(ujj,U)1)). 


The  seventh  and  eighth  terms  of  Eq . (8)  pertain  to  the  bi -frequency 

plane  in  general.  The  seventh  term  is  the  response  at  frequency  (U>1+U)2)  » 

that  is,  G^tu^jU)  ) expresses  the  normalized  response  in  the  sum  frequency 

due  to  non-linear  interactions.  Similarly,  the  eighth  term  involves 

response  at  frequency  (ui  -tu  );  that  is,  G (ut,-iu  ) is  the  normalized 

* i s 1 s 

response  in  the  difference  frequency.  In  terms  of  Figure  la, the  octant  of 
the  bi-frequency  plane  above  the  positive  axis  corresponds  to  the  por- 
tion of  the  quadratic  frequency  response  function  which  defines  sum- 
frequency  interactions,  and  the  octant  below  the  positive  axis  corres- 
ponds to  the  portion  of  the  function  which  defines  difference  frequency 
interact  ions . 


Equations  (7)  and  (8)  together  afford  a di  rec  t way  to  proceed 
in  the  present  investigation.  In  particular,  what  is  required  for  a 
hydromechan ica 1 solution  for  the  quadratic  frequency  response  is  to 
solve  the  hydromechanical  problem  for  dual  regular  wave  excitation. 
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HYDRODYNAMIC  CALCULATIONS 

A two-dimensional  method  for  evaluating  the  lateral  second-order  force 
will  first  be  developed  and  a stripwise  calculation  of  the  lateral  force  and 
moment  in  oblique  waves  will  be  performed.  Thus  the  formulas  of  a complete  set  of 
quadratic  frequency  response  functions  for  the  lateral  drifting  force  and 
moment  will  be  obtained. 

A.  The  Lateral  Second-Order  Force  on  a Cylinder  in  Transverse 
Monochromatic  Waves 

The  second-order  lateral  forces  on  a free  floating  cylinder  in  regular 
transverse  monochromatic  waves  will  be  evaluated  at  zero  forward  speed. 

The  assumptions  are  that  1)  the  contribution  of  the  second-order  poten- 
tial to  the  second-order  force  is  negligibly  small  when  compared  to  that  of 
the  first-order  potential,  and  2)  the  drift  oscillation  is  restrained. 

The  first-order  potential  which  will  be  used  here  consists  of  the  sum 
of  the  incidence,  diffraction  and  radiation  potentials.  As  is  well-known,  all 
the  potentials  satisfy  the  linearized  boundary  conditions  at  the  mean  position 
of  oscillations  of  the  body  and  the  free  surface  as  well.  Use  of  such  linear 
potentials  and  Taylor  expansions  makes  it  possible  to  estimate  the  second-order 
force  by  calculating  the  contribution  of  the  effect  of  the  instantaneous  dis- 
placement of  the  body  from  its  mean  position  to  the  resultant  time-dependent 
force. 

The  body-fixed  moving  frame  is  designated  by  o-yz  and  the  space-fixed 
frame  by  O-YZ;  the  Y-axis  rests  horizontally  on  the  mean  free  surface  and  the 
Z-axis  points  vertically  upward.  Both  frames  become  identical  at  the  mean 
position  of  oscillation  (Figure  lb). 

The  instantaneously  displaced  hull  configuration  is  identified  in 
reference  to  the  fixed  frame  by  using  the  body  coordinates  and  the  motions 
about  its  center  of  gravity. 

The  motions  of  the  center  of  gravity  of  the  cylinder  consist  of  sway, 

heave  and  roll.  The  roll -axis  passes  through  the  CG  as  shown  in  Figure  lb, 

the  coordinates  of  which  are  0,  z-z  where  z is  the  vertical  coordinate  of 

o o 

the  center  of  gravity. 
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The  displacements  are  designated  as 


S^(u>,t)  - Re | Me’,ait|  (9) 

wi  th 

s(m)  = s<m)  +!  s<m)  (10) 

where  m=2,3,l+  indicate  sway,  heave  and  roll  modes,  respectively. 

Using  the  above  displacements,  the  coordinates  of  a point  A(Y,Z)  on 
the  body  surface  at  an  instantaneous  deviated  position  are  calculated  by  the 
following  equation: 

Y = S(2>  + y cos  S(4)  - (z-zq)  s i n 

Z = + y sin  S(4)  + (z-zq)cos  . (II) 

Expanding  cos  and  sin  S^4^  in  (11)  up  to  the  second-order,  we 

have  the  coordinates  (Y  and  Z)  in  terms  of  the  displacements  and  the  body 
coord i nates , 

v . s<2>  . y . (2-zo)sW  - i v s(l*>a 

Z = s(3)  + (z-zq)  + y J S(4)  • (12) 

The  above  (12)  gives  the  amount  of  deviation  of  the  point  A, i .e .,  ( (Y-y),|_Z- (z-zQ) J | 
The  corresponding  deviation  of  the  velocity  potential  is  approximately  esti- 
mated by  Taylor  expansion  of  the  potential  at  the  equilibrium  position. 


cp(Y,Z,t)  = cp(y.z.t)  + (Y-y)cpy(y,z,t)  + |z-(z-zQ)  |cp2(y,z,t)  (13) 

where  the  <p(y,z,t)  is  the  linear  velocity  potential  which  has  been  used  in 
the  linear  ship  motion  theory. 

Use  of  the  above  potential  in  the  Bernoulli  pressure  equation  yields 
the  following  expression: 

p(Y,Z,t)  = -pgZ  - pcpt(y,z,t)  -p(Y-y)cpyt(x,y , t)  - 

-p|Z-(z-zQ)  !<Pzt(y.z.t)  - J p|vcp(y.z,t)  |a  . (14) 

The  integration  of  the  above  pressure  on  the  instantaneous  body  contour 
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gives  the  lateral  force  in  the  following  form; 

f = J p(Y,Z,t)(-dZ)  + J p(Y,Z,t)(-dZ) 


(15) 


where  dZ  = sinadS  ; dS  indicates  a segment  and  a the  angle  of  the  seg- 
ment with  respect  to  the  Y-axis,  c is  the  mean  wetted  contour,  whereas  c^’^ 
is  the  deviated  instantaneous  wetted  contour  from  the  mean  waterline. 

Use  of  Eqs.(l2)  and  (14)  in  Eq.(15)  yields  the  following  second-order 
lateral  force: 


pg8S 


(3)c(4) 


s(4)  ♦ sW  J*t 


dy 


+pJ  {ls(2)  -(z-z0)s(4)  ]'Pyt+[s(3)  + y s(4)]v2t+^  (Vcp)2}  dZ 


+pj  (gz  + cpt)dz 


(16) 


In  the  above,  B indicates  the  beam  of  the  body,  dy  = coscrds  and  dz  = s i nods , 
where  ds  indicates  a segment  on  the  body  contour  and  a the  angle  of  the 
segment  with  the  y-axis. 

The  sum  of  the  first  two  terms  is  identical  to  uu2 m where  m 

o o 

is  the  mass  of  the  cylinder.  This  is  due  to  the  heaving  equation  of  the 
cyl inder , i .e. , 

-^os(3)  = -pgBS^  - p J <PtdY  • (17) 

c 

where  the  resultant  velocity  potential,  cp  , consists  of  the  radiation 
potentials  due  to  sway,  heave,  and  roll,  and  the  incident  and  diffraction 
potentials.  However,  the  contribution  to  the  integral  in  Eq.(17)  comes 
from  the  heave-mode  radiation  and  diffraction  potentials  due  to  the  sym- 
metry of  the  body  about  its  vertical  axis. 

The  contour  Cj  of  the  last  integral  in  Eq.(l6)  consists  of  the 

wetted  contours  on  the  right  and  left  sides  of  the  hull  calculated  from 

the  mean  waterline.  Since  Z + — is  the  relative  displacement  of  the 

g t 

hull  to  the  wave  elevation,  introducing  a variable  Zp  for  the  relative 
wave  elevation  to  the  hull  displacement,  we  have  the  last  integral  in  the 
form 

P J (gZ+<?t)dZ  = -pg  J ZrdZr 

Cl  c2 

10 


mu  - m 


K-2061 


where  ca  is  tire  wetted  contour  for  the  relative  displacement  on  the  right 
and  left  sides  of  the  hull,  and  dZf  = dZ  = sinodS.  Hence, 

-pg  J zrdzr  = - j P9(r;-0  (18) 

where  r±=e,.-(S^±B/2  S^)  (19) 

[the  relative  wave  elevation  at  the  right  and  left  sides  of  the  hull] 

ex=  - j <Pt(*B/2.  0,  t)  (20) 

[the  wave  elevation  at  the  intersections  of  the  hull  and  the 
mean  free  surface] 

Substitution  of  E q s . ( 1 7)  and  (18)  into  (16)  yields  the  second-order  lateral 
force  per  unit  length  of  the  cylinder  in  the  following  form; 


f = ^moS(3)S{4)  + p J {[s(2)-(2-zo)S(4)]^yt+[s(3)+yS(4)]cpzJd 


+ \ P J (vcp)2dz  + j Pg(-r+  + rf) 


(21) 


The  above  Eq.(21)  indicates  that  the  second-order  lateral  force 
consists  of  four  components.  The  first  is  identical  to  the  gyroscopic 
coupling  term  between  roll  and  heave  motions  in  the  Eu ler-heave-equat ion. 

The  second  is  equa  1 to  the  sum  of  two  component  forces,  i.e.,  the  rate 
of  change  of  the  resultant  lateral  hydrodynamic  force  per  unit  lateral  dis- 
placement times  the  resultant  lateral  displacement  and  the  rate  of  change  of 
the  resultant  lateral  hydrodynamic  force  per  unit  vertical  displacement 
times  the  resultant  vertical  displacement.  The  third  is  the  force  due  to 
Bernoulli's  quadratic  pressure.  The  fourth  is  due  to  the  relative  wave 
elevation  on  the  hull  sidewalls. 

The  time  average  of  f in  Eq  - (2 1 ) is  expected  to  be  nearly  identical  to 
Maruo's  formula 

f " 7 pga2 Ia_ I2  (21a) 

where  a is  the  incident  wave  amplitude  and  A_  is  the  vector  sum  of  the 
wave  amplitude  ratios  due  to  radiation  and  diffraction  in  the  far  field 
Y — * -00  . 

The  present  formula  has  been  derived  by  direct  integration  of  the  pres- 
sure distribution  on  the  hull  surface  whereas  Maruo's  Eq.(21a)  is  due  to 
the  integration  of  the  pressure  in  the  far  field.  The  formula  (21),  the 
near-field  solution,  has  an  advantage  for  a clear  analysis  of  the  second-order 
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force,  though  it  does  appear  to  be  complicated.  It  serves  also  for  the 
evaluation  of  the  quadratic  frequency  response  function  for  the  lateral 
drifting  force  and  moment. 


B . The  Lateral  Second-Order  Force  and  Moment  of  a Ship  in  Oblique 
Monochromatic  Waves 

The  second-order  lateral  force  and  moment  of  a ship  at  zero  speed  in 
oblique  monochromatic  waves  are  evaluated  by  a stripwise  integration  of  the 
section  force  f in  Eq.(2l). 

The  oblique  wave  h is  given  in  the  body  coordinate  system  in  the 
fol low i ng  form 

h = a cos(v  x cosn  + vy  sinp  - uut)  (22) 

where 

a = wave  amplitude 
v = uu2/g  = wave  number 


H = angle  of  incidence  (zero  in  the  following  wave  and  n/2  in  the  beam 
wave 

/ \ 

The  linear  response  motions  ’ of  the  ship  are  determined  by  solving 
the  coupled  heaving-pi tching  and  swaying-roll ing-yawing  equations.  They  ore 
denoted  by 


. ( m)  _ - ( m)  , . - i uut 

>'  ' = Re  S'  (u)e 


s (m)  = s(m)  + iSH 


(23) 


with  m=2,3,4,5,6  (=  sway,  heave,  roll,  pitch,  and  yaw,  respectively). 

Then  Eq .(21)  gives  the  lateral  second-order  force  f(x)  per  unit  length 
of  a section  at  x“x  in  the  following  form 


f(x)  * uuamQ(x)  S^(x,0))  S^(w) 

+ P j )[s(2)(x,uu)  - (z-zq)  S(4)(u,)]cpyt(uu)dz 

+ p J j_S^(x,w)  + yS^4^  (ui)  J(f2t(a))dz  + £■  j*  I vcp(a))j2dz 
c(x)  c(x) 

+ \ P9  j_-r2  (x,u>)  + ra(x,u>) J (24) 
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wi  th 

S^(x,w)  - $^(u)  - x (in) 

S(2)(x,u))  - S*2*(u>)  + x S(6*(u>)  . 

Integration  of  the  above  sectional  force  over  the  hull  length  L yields 
the  resultant  force  and  moment  in  the  following  non-dimensional  forms; 


C,  - 


I 


F ~ I _ _3“  l f(x>dx 


T P9« 


7 PgL  L 


* Jdx  £M  { ( Jdx  J(l-«  , p itfl  dp 

L i ex 


L c(x) 


(3) 


c(x) 


a V « ' a 
L c(x) 


dz 


, » „ r.  r rvvMl*  I r I r+(x-a,)f  r-(x  ,<a).2 
2 P [dXc{x)L  a J d*  + 2 M Jdx  [-(-4-J  <y-—)  j}  (26) 


:M  = 1 “7  JxF(x)  dx 

2 Pga^L  L 


I 


^^-dx 


[ P9L3 


(3) 

a 


♦ Jdxx^i  J J(l.,o 

L C(x)  L c(x) 


V H 

)p  _X- j2 


JdxxiHl^  J p^dP.si^i,  , j yp!^Jp 

L c(*)  L C(x) 

f’f"  c{  Ir^J d*  * i P9  {*  H—V]}  . 


(27) 
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C . The  Quadratic  Frequency  Response  Function  for  the  Lateral  Drifting 
Force  arid  Moment  in  Dual  Oblique  Waves 

The  dual  wave  is  defined  by  the  sum  of  two  monochromatic  incident 
waves  (Eq.22) 

h = H ancos[v  x cos  pi  + vny  sing,  - u)nt]  . (28) 

n=l  ,2 

The  corresponding  dual  response  motion  and  potential  are  also  defined  by  the 
sums  of  the  response  to  each  monochromatic  wave,  such  as 

£ S(m^(wn)  and  D cp  (w  ) . (28a) 

n-1,2  n n=l ,2  1 " 


Substitution  of  the  above  dual  harmonic  into  Eqs.(26)  and  (27^  yields 
the  quadratic  frequency  response  functions  for  the  lateral  drifting  force  and 
moment  in  dual  oblique  waves  in  the  following  non-dimensional  form: 


CF  = 


2 P9L 


S^(w  ) 
x nr 


m K 1 , 2 


mo(x) 


n=Y,  2 


S(3)(x,oun) 


dx 


+ [dx  D 
L m=l,2 


(x,u)  ) 
' ’ nr 


m 


«{ 


x) 


n=l , 2 


<p  . (u)  ) 
yt v n7 


dz 


S^C)  CD  (uu  ) 

S — J-dx  J-  (2-z0)  s 

m=l , 2 m L c(x)  n=l ,2  an 


*{ 


dx  S 


S<3>(x.»)  cp  (lD  ) 

' nr  p £ vztv  n7 


+ Li 


m=l , 2 am  c (x)  n=l , 2 an 

c(^)  /(JO  ) „ cp  ) 

V a v m)  p.,„  P ..  v>  Yztv  n7 


dz 


/dx  J y ^ 


dz 


m=I,2  L c(x)  n=l ,2  dn 


+ 1 !d*  J ^ 


ycf(uu  ) Vcp(cu  ) 

m £ ' n7 


2 L c“(x)  m=l , 2 am  n=l,2  an 


dz 


„ r r (uu  ) r (u)  ) r (u)  ) r (cd  ) , ) 

+ £f/dx{-  x — — — s -L-2-+  ^ jljzL  s -i  , rl 

L m=l , 2 am  n=l,2  a°  m=l , 2 am  n=l ,2  an  J 


(29) 
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1 / y'  m i'  3 


S(3)  (x,iu  ) 

^ mn(x)x  2 dx 


'M  I ,a  I ^ , a *J  n "o ' ' , . a 

2 PgL  ' m»l ,2  m l n=*l ,2  n 


(x,w  ) cp  (uu  ) 

+ Jdx  x i!  _ 2-  / p s -Xl-i-  3* 


L m=l,2  m c(x)  n=1 ,2 


£ 


s<4V) 

N m 


- r v,  <Pvt(<0 

f dx  x J (z-zQ)  p £ -E dz 


m=l,2  m L c(x)  n=l,2  n 


s(3)(x,«)J 


+ dx  x E 

L m=l,2  am  c(x)  n=l,2  an 


‘ CD  (u)  ) 

r V'  Yztv  n7 

J ^ — r dz 


+ Zj 


S<4>(u,) 

m fdx 


„ <P  (ll>  ) 

f _ V1  YZtV  n7 


x J y p 


m=l,2  am  (_  c(x)  n=l,2  3n 


dz 


2 P Jdx  x J E 


vtP(U)m) 


ml  v 


v<p(u>n) 


dz 


L c(x)  m=I,2  m n=l,2  n 


I-  r*.  f ' ♦ 'nr  v n 


+ jPg  /x  x Zj 


,(uu  ) r (uu  ) r (ou  ) r (uu  ).  ) 

+ m7  v ~v  n7  . ^ jnl  p + ' n7  \ \ 


m-1,2  m n=l,2  n m= 1 , 2 m n=l,2 


(30) 


Equation  (29)  or  (30)  consists  of  seven  terms,  each  of  which  with  omission 
of  i pgL  or  j PgL*2  is  expressed  in  the  following  form; 


'k  = Rel  % \^e  '*mt}Re{  ^ Sk<"n>e~' 


-i^nt\ 


m=l,2  ' ' n=l , 2 

where  A^.B^  = complex  frequency  response  amplitudes. 

For  the  analysis  of  the  quadratic  frequency  response  function,  one  de- 


(31) 


fines  a fundamental  function  H.  . for  each  | in  the  form 

I i K K 


H.  . (uu  ,uu  ) = A.  (uu  ) B.  (uu  ) 

1 ,k'  m’  n7  k'  m7  kv  n7 


H,"  , (u>  , -u)  ) = H,  , (-U)  ,uu  ) (32^ 

1 ,k'  m’  n7  1 ,kv  m’  n7 

where  the  asterisk  indicates  the  complex  conjugate  of  the  function. 


L. 


h 

M 


► A 


. .x  MMi . 
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Use  of  the  above  functions  for  theklt'  term  of  Eqs.(29)  or  (30)  yields 
the  following  expression; 


T |Re  lHl,k(,VaJ  + Hl|kK  •"%»>] 
+ Re  Kk(v%,)«"i2lUmt] 


Re  LHi>kK'u)n)e’,2u,ntJ 


+ RC  LH|.k(u,m*UJn)  + H|,k^n.u)m)je 

V . -I  - ' c 

+ Re  LHl,k(“’in’"uVi)  + H,;k(a'n*-%1)Je 


Now  one  defines  the  quadratic  frequency  response  function  G (uu  uu  ) • 

2 N m n 

,7  . 

G2(u'm*u,n)  " 2 K .k^1^  + H1 ,kK*“tn>  J 

,7 

G2^m’'a’n^=  2 + H| , k^n» '“m)  J • (34) 

Use  of  the  definitions  (Eq.32)  makes  the  G2(w  u>  ) in  Eqs . (34)  satisfy  the 
symmetry  relationships  previously  indicated  in  Eqs. (5)  and  (6). 

Thus  the  non-dimensional  drifting  force  CF  (Eq.29)  and  moment  (Eq-30)  are 
expressed  in  terms  of  the  quadratic  frequency  response  function  G2  (Eq.?^)  in 
the  following  form; 


2PtyL  j 1 {G2<W  + 


+ 7Re  {G2^m^m>e‘i2aJmt} 


+ 2 Re  lG2(  Vn>e 


_ 1 2u.nt  ! 
J 


+ Re  k,(u»  ,vi  )e"i(UJm+a,n)tl 
*■  *■  m n J 

+ Re  |g2  -uj  )e_l  ^u,m-tl,n).t|  ( 


where  the  G2  represents  the  quadratic  frequency  response  function  for  the 
lateral  drifting  force  or  moment. 
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In  dimensional  form,  they  are  giver  by 


Mj  Ham3G2<VH*,B>  4 anG2^n‘-UJn)} 

+ j Re  {a^  G2(%1.^(n)e''2amt} 


l „ fa 


+ j Re  >a  " G (w  )e 
2 t n 2 n n 


-i2ujnt 


+ Re  {amanG2(VU,n)e",(a‘mfan)t} 

+ Re  {amanG2(V-^)e‘,(Xnran)tj 


(36) 


The  mean  drifting  force  and  moment  deduced  from  Eg. (35)  are  given  by 


CF 

J 

1 

PgL 

1 

i CM 

I 

G0(u  ,-u  ) 
2 4 nr  nr 


pgl‘ 


(37) 


where  suffixes  F and  M indicate  the  G„  for  force  and  moment,  respectively 


D.  Calculation  of  the  H,  , Function 

1 , k 


The  ^ functions  which  have  been  defined  in  Eq.(32)  are  evaluated 
in  Appendices  A,  8,  C,  and  D. 

The  Hj  ^ functions  for  the  drifting  force  are  shown  in  Appendix  A 
whereas  those  for  the  drifting  moment  are  in  Appendix  B. 

The  resultant  velocity  potential,  hydrodynamic  pressure  and  velocity 
components  are  evaluated  in  Appendix  C. 

Formulas  in  Appendices  A and  B contain  the  terms  of  hydrodynamic  pres- 
sures and  velocity  components.  These  are  evaluated  in  Appendix  C. 
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The  resultant  velocity  potential  (Eq.C-15)  is  utilized  in  calculating 
the  hydrodynamic  pressure  ( Eq - C - 1 6 to  Eq . C - 1 8)  and  the  velocity  components 
( Eq . C ~ 1 9) • It  is  also  used  in  the  evaluation  of  the  diffracted  and  radiated 
wave  amplitudes  in  the  far-field  which  serve  for  the  evaluation  of  Maruo's 
formula  (21a).  The  asymptotic  expression  of  the  potential  in  the  far  field 
is  given  in  the  Appendix  of  Reference  11. 

The  velocity  components  due  to  radiation  and  diffraction  are  calculated 
by  using  the  formulas  for  the  unit  source- i nduced  velocity  components,  which 
are  given  in  Appendix  D. 
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NUMERICAL  CALCULATIONS 


A . Comparison  of  Analytical  and  Experimental  Estimates  of  the  Mean  Drifting 
Force  and  Moment 


The  previously  defined  non-dimensional  mean  drifting  force  and  moment 

coefficients  Cr  + £ (Eq.37)  are  the  values  of  the  quadratic  frequency  re- 

r M 

sponse  functions  along  the  line  , and  it  is  this  particular  section 

of  the  functions  for  which  there  are  a number  of  conf inning  data. 


The  models  chosen  for  the  comparisons  are  1)  the  cylinder  tested  by 

hot erayama , ' ^ and  2)  the  Series  60,  0.60  C parent  which  was  used  by  Lai; 

B 

The  principal  particulars  of  the  models  are  given  in  Tables  1 and  2. 


The  non-dimensional  coefficients  Cr  and  £ (Eq.37)  are  expressed  in  terms 

F M 

of  the  previously  defined  Hj  ^ functions  (Appendices  A and  B) , each  of  which 
consists  of  four  components  as  shown  below; 


wi  th 


n = l , 4 

V 


'Fn 


n-l  ,4 

Mn 

1 

Re 

8 , 

H,  j (oj,  -uo) 

pgL 

j _ 

D 0 

M 

H,  | , (u),-uj) 

pgL‘ 

r\e 

I A F 

PgL  Rc^2 


PgL  1-2 


Re_£  H*  x(uj,-uu) 


(38) 


(39  a) 


(3  9 b) 
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-V  Re  M1.  ,(U),-U>) 

pgL  l.o 


: Re  l/'  Aw,-*') 

M3  pgLJ  ll6 


(39c) 


Cc  “ Re  H1  7(uj,-u)) 
F(|  pgL  I./ 

C — ■ “ Rc  h”  7K-w) 

M4  pgL3  ' ’ 7 


(39d) 


where  superscripts  F and  M to  the  M,  . indicate  the  functions  for  drifting 

I , k 

force  and  moment,  respectively. 


The  first  coefficients  (Eqs . 39a.  A-1*,  B-4)  arc  known  as  the  gyroscopic 

21 

coupling  terms  in  Euler's  equation  for  a body  moving  in  a vacuum.  Vossci s 

qives  a discussion  of  the  C„  describing  its  influence  on  the  steering  of  a 

Mi 

ship. 

Equation  (39b)  is  the  sum  of  two,  i.e.,  a)  the  product  ol  [the  rate  ot 
change  of  the  resultant  lateral  hydrodynamic  torce  (moment).]  and  [the  result- 
ant lateral  displacement J , and  b)  the  product  o(  [the  rate  of  change  ol  the 
resultant  vertical  hydrodynamic  force  (moment)]  and  [the  resultant  vertical 
d i sp I acement  I . 


Tlie  third  terms  are  due  to  the  Bernoulli  quadratic  hydrodynamic  pres- 
sure. It  is  to  be  noted  that  Bernoulli's  mean  second-order  torce  and  moment 

appear  as  one  of  the  four  components  of  the  total  drifting  force  or  moment. 

The  fourth  coefficients  are  due  to  the  relative  wave  elevation  along 
the  intersection  of  the  mean  free  surface  and  hull  surface  at  its  mean 
pos i t i on . 

According  to  the  aforementioned  remarks,  the  four  components  may  lie 

designated  as:  1)  the  gyroscopic  coupling,  2)  the  lorce  (moment)  rate  x dis- 

placement, 3)  the  Bernoulli  quadratic,  and  4)  the  relative  wave  components. 
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Cy  I i rider  ; 

The  four  components  Cc  of  the  mean  drifting  force  coefficients  on 

17  h"’ 

the  cylinder  versus  the  non-dimensional  wave  frequency  parameter  vb/2  are 
shown  in  Figure  2a.  Two  arrows  pointing  downward  indicate  the  resonance 
frequencies  of  roll  and  heave.  All  the  components  indicate  peaks  at  the 
roll  resonance  frequency.  The  second  peaks  ol  Hie  components  correspond!' nq 


to  the  heave  resonance  fall  slightly  apart  from  the  resonance  frequency. 
However,  the  peaks  of  the  gyroscopic-coupling  and  the  force  rate  x displace- 
ment terms  fall  nearer  to  the  resonance  frequency  than  the  others.  This  is 
attributed  to  the  fact  that  the  former  terms  arc  more  sensitive  to  the 
response  motions  than  the  latter  . The  influence  of  the  resonance  motions 
on  the  relative  wave  component  is  illustrated  in  Figure  2b. 


The  gyroscopic-coupling  term  C is  negative  all  along  the  frequency 

f l 

axis.  This  is  due  to  the  fact  that  for  the  model  the  phase  differences  be- 
tween the  rolling  and  heaving  displacements  are  greater  than  n/2  in  the 
entire  frequency  domain. 


The  force  rate  x displacement  term  indicates  its  peak  at  roll  resonance 
and  negative  peak  in  response  to  heave  resonance. 


The  Bernoulli  quadratic  component  appears  to  be  negative  in  the  entire 
frequency  domain,  which  results  from  the  force  unbalance  due  to  the  quadratic 
pressure  asymmetric  with  respect  to  the  vertical  z-axis. 

The  relative  wave  term  Cr  is  dominant  among  the  components.  It  exceeds 
and  counteracts  all  the  rest  of  the  components  and  thus  yields  the  positive 
resultant  drifting  force,  which  remains  below  unity.  This  confirms  the 
calculation  and  agrees  with  the  energy  conservation  law.  The  behavior  of  the 
magnitudes  of  the  ratio  of  relative  elevation  to  Incident  wave  amplitude 
|r(/a|  is  illustrated  in  Figure  2b,  where  A indicate  the  locations,  i.c., 
the  intersections  of  the  hull  surfaces  at  equilibrium  with  the  mean  free 
surfaces  along  the  positive  and  negative  y-axis.  When  the  body  is  free,  it 
is  evidently  affected  by  the  resonance  motions. 

The  asymptotic  behavior  of  the  force  components  lor  the  infinite  fre- 
quency are  estimated  as  follows: 


As  the  frequency  approaches  infinity,  the  response  motion  vanishes  and  it 
results  that  I)  the  hull  motion  dependent  gyroscopic-coupling  and  force  rate 
x displacement  terms  will  vanish,  and  2)  the  relative  wave  elevations  for 
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free  body  reach  to  the  state  of  perfect  reflection  against  a fixed  body. 
In  mathematical  form  we  have 


f requency  oo 


(40a) 


and 


|r_|  - 2a 


(40b) 


Hence,  referring  to  ( 39d)  and  (A-28) , 


According  to  the  energy  conservation  in  perfect  reflection  at  infinite 
frequency,  we  have 


Thus,  it  is  inferred  that 

<L  - -1.0 

■"a 


(40c) 


(40d) 


(40c) 


For  the  vanishingly  low  frequency,  the  four  components  are  estimated  in  the 
fol lowi ng  manner  : 


The  phase  difference  between  heave  and  roll  is  n/2  in  Eq.(A-4),  thus 

Cp  -•  0;  since  the  hydrodynamic  force  rate  is  linearly  proportional  to  the 

wave  frequency  (Eq.A-6) , Cr  -•  0.  Since  the  heave  phase  is  zero  and  roll 

'"a  (3)  (M 

phase  is  n/2,  the  components  S'  ' and  S'  vanish  together  with  the  terms 


s c 

in  Eq.(A-30).  Thus  due  to  (A-28),  C 


containing 
mission  of  energy  is  perfect,  C 


0.  Since  the  trans- 


0.  Hence  it  is  inferred  that  the  Bernoulli 


quadratic  term  C 


F3 


0. 


Figure  2c  illustrates  a comparison  of  the  analytical  methods  and  experi- 
mental values.  Maruo's  far-field  theory  is  presented  by  a dotted  line  and 
the  present  near-field  method  is  indicated  by  a solid  line.  Both  results  arc 
in  good  agreement  as  was  expected.  The  near-field  value  indicates  the  effect 
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of  roll  resonance,  whereas  the  far-field  procedure  does  not.  The  experi- 
mental values  indicate  higher  values  than  unity  in  the  high  frequency  region. 


Ship; 

The  four  constituents  of  the  beam  sea  mean  drifting  force  acting  on 
the  ship.  Reference  12,  are  illustrated  in  Figure  3a  in  coefficient  form. 

The  behavior  of  the  following  coefficients  appear  to  be  similar  to  that 
of  the  cylinder  in  that  1)  the  hull -mot  ion  dependent  terms  C_  and  C_  show 

•"1  r2 

their  peaks  at  (or  near)  the  heave  resonance  frequency  (see  Figure  3b); 

2)  the  phase  difference  between  heave  and  roll  as  shown  in  Figure  3b  appears 
to  be  greater  than  tt/2,  resulting  in  a negative  gyroscopic-coupling  term; 
and  3)  the  relative  wave  and  Bernoulli's  quadratic  terms  reach  to  their 
peaks  in  the  neighborhood  of  the  heave  resonance. 

The  resultant  mean  drifting  forces  and  moments  are  illustrated  in 
Figure  3c.  Maruo's  method  and  experimental  model  data  are  compared  with 
the  results  of  the  present  calculation.  The  comparison  shows  excellent 
agreement  in  beam  seas.  In  oblique  seas,  the  two  analytical  methods  are  in 
good  agreement,  as  shown  in  Figures  3d  and  3e. 


The  Irregular  Frequency  Phenomenon: 

The  presence  of  a mathematical  discontinuity  in  the  mean  drifting  force 

at  vB/2  «s  1.87  for  the  cylinder  has  been  investigated.  The  discontinuity 

is  due  to  the  ei genf requency  of  the  fluid  motion  inside  the  hull.  Referring 

1 8 

to  the  work  of  Ogilvie  and  Shin,  a modification  of  the  Green's  function  has 
been  performed  just  for  the  heave  mode  since  its  contribution  is  regarded 
dominant.  After  a trial -and-error  study,  we  have  selected  the  constant  C 
of  0.033. 

Figure  4 illustrates  that  the  modification  is  virtually  effective  in 
reducing  the  discontinuity  behavior  of  the  drifting  force. 


23 


R-206 1 


B . The  Quadratic  Frequency  Response  Function  for  Lateral  Force 

The  purpose  of  the  development  of  computational  methods  for  quadratic 

response  functions  is  to  enable  improvement  in  the  capability  of  predicting 

the  non-linear  lateral  force  and  moments  induced  by  i rregular  waves.  In 

practical  terms,  the  important  part  of  the  sway/yaw/surge  response  of  ships 

is  not  the  part  which  may  be  predicted  by  linear  methods,  but  the  low- 

frequency  non-linear  part;  that  is,  the  forces  which  occur  at  frequencies 

below  the  lowest  wave  excitation  frequencies  of  significance.  Given  the 

quadratic  frequency  response  function  and  a specification  of  the  irregular 

waves,  the  contribution  to  lateral  forces  of  the  quadratic  non-linearities 

2 

treated  herein  may  be  computed  either  in  the  time  or  frequency  domain. 

Insofar  as  engineering  predictions  are  concerned,  it  appears  reasonable  to 
consider  the  predicted  very  low  frequency  forces  as  excitation  to  a linearized 
representation  of  sway  and  yaw  (as  in  Reference  3)  so  that  analytical  predic- 
tion of  the  important  features  of  the  non-restored  modes  of  motion  in  irregular 
seas  is  within  reach  - to  the  extent  of  course  that  the  estimated  quadratic 
frequency  response  functions  are  valid. 

The  portions  of  the  quadratic  frequency  response  function  pertaining  to 
interactions  between  two  different  frequencies  are  of  importance  in  the  above 
respect.  Since  the  preceding  sections  of  the  report  have  indicated  some  quite 
favorable  exper imenta l/ana lyt i cal  comparisons  for  the  mean  drift  and  moment 
function  (Ga  (w,-u>))  , the  next  logical  step  would  be  to  compare  analysis  and 
experiment  for  the  general  function  (Ga  (“1  ,wa) ) , (or,  at  least  parts  of  it). 
Because  no  suitable  experiments  had  been  conducted  as  of  the  time  of  writing, 
this  next  logical  step  was  not  possible.  Thus  the  most  that  was  reasonable 
to  do  in  this  circumstance  was  to  compute  a typical  quadratic  response  function 
for  assumed  experimental  conditions,  and  to  examine  its  plausability  relative 
to  previous  experience. 

To  this  end  the  model  ship  (Series  60,  C =0.6)  as  used  in  Reference  12, 
and  for  which  the  mean  drift  force  response  was  considered  earlier,  was 
assumed  to  be  at  a sixty-degree  heading  to  waves,  at  zero  speed,  and  the  com- 
putations indicated  in  Eqs.(29,  (31),  (32),  and  (34)  were  carried  out  for  the 
lateral  drift  force  components.  The  values  of  the  function  were  non-dimen- 
sionalized  in  the  form  implied  by  Eq.(37);  that  is,  the  dimensional  values  of 
(forcej/length)"  were  divided  by  pgL.  It  was  convenient  to  noa-dimens iona I i ze 
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the  frequencies  in  the  following  form; 


where : 


o 

m 


U)  /id 

n IL 


( = v'L/X  ) 


■ V2ng/L 

L = ship  length 
X = wave  length 


The  computation  was  carried  out  for  all  combinations  of  a non-dimensional 
input  frequency  range  from  0=0.373  to  2.2  (X/L  = 0.2  to  7.2).  Thirty 
equally  spaced  frequency  steps  were  utilized  within  this  range  to  reasonably 
define  the  function. 

In  order  to  visualize  the  results  it  is  advantageous  to  utilize  the 
planes  of  symmetry  discussed  in  conjunction  with  Figure  la, and  thus  also  of 
advantage  to  introduce  the  frequency  transformation  of  Reference  2.  This 
transformation  maps  the  ,u^  plane  (see  Figurela)  into  a sum  and  difference 
frequency  plane  in  accordance  with 

-uug  = difference  frequency 

Q3  = u^+uug  = sum  frequency  (41) 

The  0^  ,Qa  plane  is  shown  mapped  into  the  .uig  plane  in  Figure  5.  The 

,Qa  axes  are  coincident  with  the  lines  of  symmetry  noted  in  the  discussion 
of  Figure  la.  In  accordance  with  the  discussion  of  Figure  la,  the  fundamental 
symmetry  of  the  quadratic  frequency  response  function  allows  complete  defini- 
tion if  the  function  is  known  in  a quadrant  of  the  plane  defined  by  the  lines 
of  symmetry.  The  quadrant  of  the  QltQs  plane  (Figure  5)  in  which  both  sum 
and  difference  frequencies  are  positive  corresponds  to  the  quadrant  selected 
earlier  in  interpreting  the  response  to  dual  harmonic  excitation.  Also  in 
accordance  with  the  earlier  discussion  of  dual  harmonic  excitation,  the  sum 
frequency,  Qa,  corresponds  and  is  numerically  equal  to,  the  output  frequency. 
Any  combination  of  input  frequencies  and  uvj  which  are  on  the  line  n3  = con- 
stant will  generate  a component  at  that  frequency.  For  example,  the  mean 
drift  force  operator  involves  output  frequency  of  zero  and  is  located  on  the 
Ql  axis.  The  frequency  transformation  involves  no  distortion  of  the  magnitude 
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of  the  quadratic  frequency  response  function,  merely  providing  a somewhat 
more  convenient  way  of  looking  at  it.  Since  the  non-dimensional  frequencies 
(on)  differ  from  the  dimensional  frequencies  (u^)  only  by  a scale  factor, 
and  Oa  may  be  substituted  for  and  in  Eq.(^l)  and  the  resulting  ,fl2 
system  taken  to  be  non-dimensional. 

Figure  6 is  an  isometric  plot  of  the  modulus  of  the  analytically  de- 
termined quadratic  frequency  response  function.  The  "view"  is  toward  the 
origin  of  the  bi-frequency  plane  from  a position  above  the  axis  (Figure  5). 
(The  difference  frequency  axis  (fy)  is  to  the  left,  the  sum  frequency  axis  is 
to  the  right,  and  the  modulus  of  the  function  is  "vertical".)  The  lines  de- 
fining the  function  contours  are  in  the  nature  of  section  lines  which  are 
parallel  to  the  and  Q2  axes.  Each  intersection  corresponds  to  an  analyti- 
cally computed  point,  and  the  connecting  line  segments  are  straight.  There 
were  actually  about  twice  as  many  points  computed  as  are  shown.  In  the 
figure  the  plotted  domain  of  the  bi-frequency  plane  appears  triangular.  This 
is  the  result  of  computing  the  function  to  a maximum  absolute  input  frequency 
of  2.2.  The  right-to-left  section  at  the  front  of  the  plot  is  indicated  with 
vertical  lines  defining  the  position  of  computed  points.  This  section  repre- 
sents the  modulus  of  the  function  along  the  line  = 2.2.  There  is  an 
exactly  flat  portion  shown  in  the  center  of  the  surface  because,  as  an  economy 
measure,  no  computations  were  carried  out  for  input  frequencies  between  0 and 
0.373,  and  for  plotting  purposes  the  magnitude  of  the  interactions  between 
these  frequencies  and  all  others  was  assumed  to  be  zero. 

Figure  7 has  been  prepared  as  an  aid  to  discussion,  as  well  as  to  round 
out  the  presentation  in  Figure  6.  At  the  bottom  of  this  figure  the  linear 
roll  and  heave  responses  are  shown  in  non-dimensional  form  plotted  on  input 
frequency,  o.  The  lower  two  frames  of  the  figure  indicate  the  modulus  of  the 
quadratic  frequency  response  function  on  the  lines  of  symmetry;  that  is,  along 
the  lines  Q1=0  and  Cl8  = 0 . In  these  special  cases  there  is  a one-to-one 
correspondence  between  input  frequency  and  the  sum  and  difference  frequencies 
Qj,  and  Qa  since  the  function  may  in  these  cases  be  considered  to  indicate  the 
mean  shift  and  second  harmonic  response  for  excitation  at  a single  frequency. 

Qual i tat i vely,  Figures  6 and  7 show  that  the  modulus  of  the  computed 
quadratic  frequency  response  function  resembles  that  which  is  expected  from 
previous  work  with  added  resistance  (Reference  1).  There  are  two  main  humps 


shown  within  the  range  of  computation,  one  associated  with  difference  fre- 
quency interactions,  one  associated  with  sum-frequency  interactions,  and 
the  two  humps  are  not  symmetrical. 

Several  more  detailed  observations  can  be  made.  As  was  noted  in  the 
presentation  of  Figure  6,  the  apparent  flat  valley  in  the  middle  comes  about 
because  no  computations  were  made.  The  edge  of  this  valley  has  some  minor 
"cliffs"  which  (translated  into  the  o3  plane)  occur  along  the  lines 
Oa  — t 0.4  for  greater  than  about  1.  This  means  simply  that  the  inter- 
actions of  long  low  frequency  waves  with  much  shorter  high  frequency  waves 
are  a good  deal  larger  than  the  value  of  the  quadratic  frequency  response 
function  in  the  planes  of  symmetry  would  lead  one  to  expect.  Comparing 
Figures  6 and  7,  the  most  important  magnitudes  of  the  function  occur  near 
heave  resonance  so  that  it  appears  that  non-linear  interactions  between  very 
long  waves  and  waves  causing  heave  resonance  can  be  non-negl i gi ble,  though 
of  secondary  magnitude. 

The  next  observation  from  Figure  6 is  that  there  appears  to  be  two  chains 
of  secondary  (and  not  well  resolved)  peaks  in  the  modulus  of  the  function  run- 
ning parallel  to  the  "cliffs"  just  mentioned.  Comparing  the  frequencies  in 
Figure  7 it  is  clear  that  these  peaks  are  the  result  of  resonant  rolling.  As 
shown  in  Figure  7.  the  computed  rolling  is  lightly  damped  because  the  present 
computation  is  made  for  zero  ship  speed,  and  low  "viscous"  roll  damping  effects 
have  been  assumed.  Rolling  is  thus  shown  to  have  a secondary,  though  possibly 
not  entirely  negligible  effect  on  the  lateral  drift  force  function  since  the 
computations  show  interactions  between  responses  at  roll  resonance  and  at 
frequencies  approaching  heave  resonance  to  be  relatively  larger  than  the  mean 
drift  response  near  roll  resonance  shown  in  Figure  7. 

The  most  important  practical  part  of  the  function  in  Figure  6 is  that 
to  the  left,  the  part  which  influences  low  frequency  response,  and  this  part 
of  the  function  appears  smooth  and  relatively  uncomplicated.  There  is  a notch 
in  the  mean  operator,  Figure  7,  at  a frequency  of  a = 2.1  (0X  = 4.2),  and  this 
is  traceable  to  an  "irregular"  frequency.  As  noted  in  the  previous  section, 
the  technique  adopted  for  minimizing  the  effects  of  irregular  frequencies 
does  not  completely  eliminate  the  problem.  Effects  of  this  irregular  fre- 
quency are  included  in  some  way  for  all  estimates  of  the  quadratic  frequency 
response  function  involving  input  frequency  of  o = 2.1.  The  largest  input 
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frequency  for  which  computations  were  made  (o  = 2.2)  corresponds  to 
(wave  1ength)/ship  length  of  0.2.  As  far  as  showing  a picture  of  the  com- 
plete function  this  is  too  long  a wave  length,  but  any  wave  length  would 
be  too  long  since  it  was  shown  that  the  mean  drift  force  operator  is 
asymptotic  to  a constant  as  wave  length  goes  toward  zero.  It  can  only  be 
surmisec,  but  it  appears  that  the  quadratic  frequency  response  function 
continues  indefinitely  along  the  Qj  axis  about  as  shown  in  Figure  6 for  the 

highest  frequencies  computed.  This  last  feature  of  the  function  complicates, 

2 

but  does  not  prevent,  practical  predictions  for  irregular  waves  since  the 
magnitudes  of  very-very  short  waves  in  an  irregular  sea  are  ordinarily  ex- 
tremely small  relative  to  the  magnitudes  of  waves  near  ship  length  and 
longer. 

The  last  feature  of  the  computed  quadratic  frequency  response  function, 
the  most  significant  part  of  the  sum-frequency  interactions  (to  the  right  in 
Figure  6),  is  the  one  that  the  least  is  known  about.  The  computed  peak  modulus 
for  second  harmonic  response  (Figure  7)  is  about  twice  that  for  the  mean  drift 
force  response  and  has  a hump  and  hollow.  Why  these  come  about  is  not  immedi- 
ately obvious.  It  is  again  a surmise,  but  if  the  difference  frequency  part  of 
the  function  continues  indefinitely  along  the  axis,  significant  values  of 
the  sum-frequency  part  would  be  expected  to  continue  indefinitely  along  the 
Da  axis.  The  range  of  the  computation  is  too  short  to  convincingly  demonstrate 
the  behavior  of  the  function  in  this  direction. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

What  has  been  accomplished  thus  far  is  the  development  of  a "near 
field', 1 second-order  strip  theory  for  lateral  force  and  moment  on  a ship 
at  zero  speed  which  is,  moreover,  formulated  in  such  a way  that  the  results 
of  the  hydromechanical  analyses  may  be  immediately  connected  with  a long 
existing  general  input/output  theory  for  weakly  non-linear  systems.  As  a 
result  of  the  formulation  the  computation  of  all  the  quadra t ic, non- 1 inear 
interaction  responses  between  pairs  of  waves  (the  "complete"  quadratic 
frequency  response  function)  is  hardly  more  expensive  than  a careful  linear 
motions  analysis  utilizing  the  strip  "close-fit"  approach. 

To  the  extent  that  the  results  of  the  present  development  can  be 
compared  with  experiment,  the  comparision  appears  favorable.  However,  the 
experimental  data  available  are  concerned  only  with  some  very  special  values 
of  the  quadratic  frequency  response  function  (in  fact,  values  of  the  mean 
drifting  force  and  moment  measured  in  monochromatic  waves).  Thus  while  the 
present  analysis  is  capable  of  predicting  what  appears  to  be  a plausible 
complete  quadratic  frequency  response  function,  experimental  confirmation, 
either  direct  or  indirect,  of  the  entire  computation  is  lacking. 

In  order  to  make  further  progress  in  this  area  some  experiments, 
which  go  beyond  the  convent iona 1,  need  to  be  carried  out,  with  the  objectives: 
a)  of  enabling  comparisons  of  experiment  and  theory  in  at  least  the  most 
important  areas  of  the  bi-frequency  domain,  and  b)  of  enabling  comparisons 
to  be  made  between  lateral  drifting  forces  and  moments  in  irregular  seas  and 
those  predicted  completely  by  analytical  methods. 

Since  the  influence  of  the  quadratic  non-linearities  in  lateral  wave 
induced  forces  can  be  important  in  applications  where  precise  directional 
control  of  ships  is  required,  an  extent  ion  of  the  present  methods  to  include 
forward  speed  is  also  recommended. 
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TABLE  1 

PRINCIPAL  PARTICULARS 

The  particulars  of  Koyerayama's  model,  which  have  been  used  in  the  present 
calculation,  are  tabulated  below; 


Cylinder  Shape  of  a Lewis-Form  Section 


Section  Area  Coefficient 
(due  to  close-f i t) 

0.9326 

Beam,  ft 

1.64 

Draft,  ft 

0.656 

Rol 1 Gyrad ius , ft 

0.5711 

Metacentric  Height,  ft 

0.12464 

Vertical  Center  of  Gravity 
(below  water  1 i ne) 

-0.07544 

Roll  Period,  sec 

1.895 

Roll  Damping  Factor 

0.62 

original  value;  0.942 
’original  value;  -0.108 


I 
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TABLE  2 

MODEL  PARTICULARS 

Series  60,  C =0.6 
B 


LBP,  L,  ft 

5.00 

Breadth,  B,  ft 

0.667 

Draft  (level  trim),  T,  ft 

0.267 

Displacement  (F.W.),  A,  lb 

33.27 

LCG  (abaft  midsection),  ft 

0.075 

VCG  (below  water  1 i ne) ,0G,  ft 

0.022 

Rudder  area,  fts 

0.030 

Waterplane  Area,  Aw,  ft2 

2.355 

Load  Waterline  Coefficient,  C 

w 

0.706 

Pitch  Gyradius,  x^,  ft 

1.275 

Yaw  Gyradius,  x , ft 

X 

1.275 

Roll  Metacentric  Height, GM  ,ft 

m 

0.025 
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FIGURE  la.  SYMMETRY  OF  THE  QUADRATIC 

FREQUENCY  RESPONSE  FUNCTION 
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FIG.  3b.  SWAY,  HEAVE  AND  ROLL  RESPONSES  OF  THE 
SHIP  MODEL  IN  BEAM  SEAS  (yx=  90°  ) 
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MOMENT  ON  THE  SHIP  MODEL  IN  BEAM 


FIG.  3d.  THE  MEAN  DRIFTING  FORCE  ON  THE  SHIP  MODEL  IN  QUARTERING  SEAS 


MARUO'S  THEORY 
PRESENT  METHOD 


THE  MEAN  DRIFTING  FORCE  AND  MOMENT  ON  THE  SHIP  MODEL  IN 
QUARTERING  SEAS  (fi  = 30°) 


FIGURE  4.  THE  IRREGULAR  FREQUENCY  PHENOMENON  IN  THE  MEAN  LATERAL  DRIFTING 
FORCE  COEFFICIENT  OF  A CYLINDER 
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APPENDIX  A 


THE  H|,k  FUNCT,0NS  F0R  THE  LATERAL  DRIFTING  FORCE 


1 ) Hl,l  for  'l 


§(4) (m  ) 

i -*Ur.2— ' ’ 


-(3) 

f <*  (x*a,n^  - i uun  t 

e | Re  I ! dx  m (x)  uu  e n 

J LJl  °Wn=l,2  n an 


H,  . (u)  ,uu  ) = 

I , I m n7  a 


S^(wJ  „ _ S^(x,u>) 


J dx  m (x) 


m L 


(A  - 1 ) 
(A  — 2) 


Since  S^(x,u))  = S^H  - x S^(u>)  , 

J n’(x)dx  = m (ship  mass)  ; |’x  m (x)dx 

L L ° 


= 0 


(A-3) 


we  have 


„ S‘(3)  (x,uu  ) §(3)  (a)  ) 

J dx  m (x)uj  3 2-  = ou  2 m " 

¥ o'  7 n a - 


no  a 


thus , 


Re  H.  . (uu  ,iu>  ) = 
1,1  m n7 


s(4)(^) 

n-  ' rr\‘ 
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n a o 


s!4)(^j 


s m'  „ a s ' n 
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— ; » M 
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2)  H |2  for  l2 


l2  = Re 


1 Idx  V 

L m*l, 


, S ^ x * am^  -•‘Vl  f „ s-  a -iuLt\  ) 

■2~  • }ReU/vr.2p^^re  ^ 


-(2) 

S (x,uj  ) .9  (00  ) 9 (uu  V 
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s[ ^ (uu  ) 9 (a)  ) 

Re  H.  = 2L  f dxp«  f -^-^L  dz 
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S(4)(w)  , tp  (w  ) 

„ ,,  c ' nr  f . p zs  rv 

Re  H1 ,5  = 3 J dx  P*n  ,J  V 5 dz 

m L c(x)  n 


(u>  ) 

s v nr 


r r V(U,n) 
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6)  H 6 for  l6 
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The  velocity  components  of  the  water  are  evaluated  by  applying  formula 
(C-19)  in  Appendix  C. 
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7)  H),7  f°r  '? 
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Formulas  (A-30)  and  (A -3  1 ) are  evaluated  by  employing  ( C - 1 7)  in  Appendix  C. 
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APPENDIX  B 

THE  H|>k  FUNCTIONS  FOR  THE  LATERAL  DRIFTING  MOMENT 
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7)  H]j7  for  l? 
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where  evaluation  of  the  r is  carried  out  by  using  Eqs.(A-30)  and  (A-31)  in 
Appendi x A. 

The  velocity  components  of  the  water  on  the  hull  surface  as  shown  in 
(g-6)  _ (b-25)  are  evaluated  by  using  the  formula  ( C - 1 9)  in  Appendix  C. 
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APPENDIX  C 

THE  RESULTANT  POTENTIAL,  PRESSURE  AND  VELOCITY 

As  indicated  in  Appendices  A and  B,  we  will  now  evaluate  the  resultant 
velocity  potentials,  the  hydrodynamic  pressures  and  velocity  components  on 
the  hull  surface.  The  resultant  potential  cp  of  the  water  flow  due  to  the 
ship  motion  with  drift  restrained  in  oblique  regular  waves  consists  of  the 
incidence,  diffraction  and  rad i at  ion  wave  potentials.  These  potential  components 
are  evaluated  in  a stripwise  manner  described  as  follows: 

Let  the  incident  angle  of  the  wave  be  designated  by  p as  shown  in 
Figure  lc  and  let  the  wave  progress  in  the  positive  Y-direction.  Then  the  wave 
prof i le  is 

h|  = a cos(v  ysinp  + v xcosp  - out)  (C-l) 

where 

a = wave  amplitude 

v = wave  number  (uT’/g) 

to  = circular  frequency  of  the  wave. 

Now,  suppose  two  vertical  control  planes  cut  the  body  at  x and  x+dx,  and  observe 
the  wave  motion  within  the  fictitiously  confined  domain  which  is  infinifcly  ex- 
tended in  the  lower  half  domain,  i.e.,  in  y-  and  z-d i rect ions . This  domain  is 
designated  the  "strip  domain."  The  wave  equation,  Eq.(C-l)  can  be  interpreted 
in  this  domain  by  noting  that  the  term  vysinp  determines  the  wave  form  in 
the  strip  domain  and  that  the  term  vxcosp  represents  the  phase  shift  of  the 
incident  wave  at  y=0  and  x=x  relative  to  a crest  at  the  origin.  Thus  the 
potential  of  the  incident  wave,  Eq.(C-l)  is  defined  in  the  strip  domain  in  the 
form 

= -*^7  e *sin(vy  sinp  + vx  cosp  - cot ) . (C-2) 


This  potential  consists  of  even  and  odd  functions  with  respect  to  y 
are  expressed  in  the  forms: 

<f|  = ^7  cV/s  i n ( vy . s i np)  • cos  (vx  cosp  - u>t) 

C Q3  V 2 

9 1 = e cos(  vy -s  i np)  • s i n(vx  cosp  - u>t) 

Cl 


and  they 


(C-3) 


' 


R-2061 


The  odd  function  is  applied  to  represent  the  asymmetric  flow  about  the  z-axis, 
while  the  even  function  is  applied  to  the  symmetric  flow.  This  functional 
resolution  will  be  utilized  in  the  kinematical  boundary  conditions  in  the 
following  section. 


Diffracted  Wave  Potential  in  a Strip  Domain 

We  assume  the  diffraction  potential  in  the  strip  domain  in  the  form 


9, 


(m) 


(x,y,z;P,  t)=Re  (s) -G(y  ,z,11,§)ds  e'  (v*cosp  an) 


(C-4) 


c designates  the  mean  wetted  contour  of  the  strip  section  and  m 


Qp  ’/ (s)  designates  the 


where 

designates  the  mode  of  excitation  Lm=2 , 3 sway,  heave] 
unknown  complex  source  intensities  distributed  over  the  strip  surface.  These 
source  intensities  are  determined  by  satisfying  the  kinematical  boundary  con- 
dition on  the  body  surface 


% 


dcp 


(C-5) 


These  source  intensities  depend  on  the  mode  of  excitation,  the  geometry  of  the 
body  and  the  incident  wave.  The  function  G(y,z,T],§)  [Reference  20]  is  the 
two-dimensional  pulsating  source  potential  of  unit  intensity  at  the  point  (T|,£) 
in  the  lower  half  plane.  The  exponential  term  e,Vxcosp  represents  the  loca- 
tion of  the  strip  x=x  where  the  disturbance  occurs  in  response  to  the  oblique 
incident  wave  of  wave  number  v and  incidence  p. 

Since,  apd  G are  the  complex  source  intensity  and  Green's  func- 

tion, respectively,  let  them  be 

qH  = (>)  + i>) 

D Dr  VD  i 


G - Gr  - IG, 


where  i-  V-!  and  where  qI"'^  and  are  real  and  imaginary  parts  of  Q 


(m) 


'Dr 


D i 


and  Gr  and  -G.  are  real  and  imaginary  parts  of  G 
Equation  (C-h)  is  then  changed  to 


JV* 


P^Gr  + Op"i^Gj/Jds-cos(vxcosp-o)t) 
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(C-6) 


C2 


R-2061 


(m) 


This  potential  <p 
condition,  Eq.(C-5),  which  becomes 


•s  determined  by  satisfying  the  kinematical  boundary 


•Q 

3n 


’n 


A (3> 


Ocf 


on 


(C-7) 


on  the  fixed  body  surface.  This  boundary  condition  is  specifically  written 
in  the  following  forms  by  dropping  the  suffix  D,  i.e.,  for  sway  (m=2) , 


S aU)i!2)  * £ u!.2?  j<2) 
j-i  j 'J  j.i 


iN+j' 


vz; 


_U)e  [s  inn-cos  (vy.  sinn)sina;  - sin(vy.  sinp,)cosa;] 


-s  - e «'2>  i!2)  = o 

j=l  J 'J  j=l  N+j  ij 


(C-8) 


and,  for  heave-exciting  (m=3) , 


" <S>\  A})  -o 

j = l J 'J  j=l  N+J  IJ 
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j = l J 'J  j=1  N-tj  ij 
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-oue  1 [s i np • s i n ( vy . s i np.)  s i no.  + cos(vyj  s i np) cosa: ] 


(C-9) 


Here  a.  is  the  angle  made  between  i segment  and  the  positive  x-axis  and  both 
I j and  J.j  formally  represent  the  normal  derivatives  -jGrds  and  ^j-jGjds, 


and  are  called  the  "influence  coeff icients"  wi th  specific  definitions  given  in 
the  Appendix  of  Reference  ; and  Q.j  ^ = and  C.^+j=  Q-  ^ at  the  j’t^1  seg- 

ment. 


If  N sources  are  distributed  over  the  strip  surface,  we  obtain  a 
(2N  x 2N)  simultaneous  equation  system  with  2N  unknowns,  i.e.,  the  real  and 
imaginary  parts  of  the  N sources. 


Since  the  right-hand  sides  of  Eqs.(C-8)  and  (C-9)  represent  the  expression 


C3 


F. 


a"  5rT’  the  so,ut,on  Qq  gives  the  source  intensity  per  unit  amplitude  of 
the  incident  wave  and  the  insertion  of  this  in  Eq.(C-6)  then  provides 
the  diffraction  potential  per  unit  amplitude  of  the  incident  wave  cp^  ^ . 

In  reference  to  (C-6) , the  complex  diffraction  potential  cp^/a  in  the  strip 
domain  is  defined  by  its  real  and  imaginary  components  cp^/a  and  cp^/a 
in  the  following  forms: 


,(m)  *(m) 
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(C-l  I) 


The  Radiation  Potential  in  the  Strip  Domain 

From  the  solutions  of  the  coupled  heaving  and  pitching,  and  the  coupled 
swaying,  rolling  and  yawing  motions  of  a ship  in  oblique  seas,  according  to 
a procedure  for  motion  calculation,  we  obtain  the  five  responses  S ^ . The 
displacements  of  the  section  consist  of  three  degrees  of  freedom:  sway,  heave 
and  roll  which  are  induced  by  the  motion  of  the  ship,  namely, 

s(2)M  -s<2>  + *s(6) 

S(3)(x)  -S<J>  -,s(5) 


*<->(,)  .S« 


The  complex  radiation  potential  <p^m^/a  which  is  due  to  the  actual 

K 

response  motion  of  the  section  in  the  incident  wave  is  given  by  the  following 
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express  ions  : 
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(C- 1 4) 


b(m) 


The  7 in  the  above  presents  the  radiation  potential  per  unit  motion  dis- 

placement,  i.e.,  = 1 and  = 0 . 1 1 i s to  be  noted  that  the  sway 

-(21c  s -(21  -(41 

displacement  S has  to  be  replaced  by  S '4ZQS  ' when  the  rolling  axis 

passes  through  the  origin  0 which  is  defined  in  Figure  I-b.  Since  the  section 

motion  displacement 

satisfy  Eq - ( C- 1 2) . 


motion  displacement  is  induced  by  the  ship  motion,  and  should 


In  view  of  the  above,  the  resultant  velocity  potential  cp/a  for  each 
section  is  given  by  the  following  expressions: 
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A'* 


The  Hydrodynamic  Pressure 

The  complex  resultant  hydrodynamic  pressure  -E-  per  unit  wave 
is  derived  from  the  potential  and  % , as  the" fo I lowing: 


pc 

- = - pu)  _ 


T = pa,T 


Substituting  Eqs. ( C - 1 5)  into  Eq . (C- 16) , we  have 
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(C- 19)  i n Appendix  C 


The  Veloci ty  Components 


The  velocity  components  are  derived  from  the  potential  Eqs.(C-15) 
the  following  form; 


e sinp  cos[  v(ys  i np  + xcosp)  ] 


cos (vxcosp) 


u)  e sin[v(ysinp  + xcosp)] 


cos (vxcosp) 


«)  e s inps  i n[v(ys  i np  + xcosp] 


(vxcosp) j 


sin( vxcosp)  + 


•u)  e ZcosLv(ysinp  + xcosp)] 


sin( vxcosp)  + 


Formulas  (A-30)  and  (A -3  I ) are  evaluated  by  employing  ( C - 1 7)  in  Appendix  C. 
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The  velocity  components  due  to  the  radiation  and  diffraction  potentials 

i and  4 In  Eqs.(C-19)  are  calculated  by  using  the  formulas  (C-10),  (C-ll) 

R D 

and  (014)  and  the  formulas  for  the  velocity  component  due  to  the  source 
potential  per  unit  intensity  in  Appendix  0.  Also,  the  symbols  x and  y 
used  in  Appendix  0 are  identical  to  y and  z . 
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APPENDIX  D 

THE  VELOCITY  COMPONENTS 
DUE  TO  THE  UNIT  PULSATING  SOURCE,y 

The  velocity  components  due  to  the  unit  pulsating  source  (G=G  r~ ' G j ) 
along  the  cylinder  contour  are  evaluated  by  applying  Frank's  closed-fit 
method. ^ The  fluid  velocity  components  at  the  i*^  midpoint  (x.ty.)  due 
to  the  source  on  a segment  Sj  for  mode  m are 
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where 


-ik(2  C) 

CO  I 

Gr  = J [log(z  -C)  - log(z.-;)  + 2P.V.  J - dk3ds 


S. 

J 


- k 


- i k(  z +C) 


(-l)m  J [log(z.+C)  - log(z.+C)  + 2P.V.  J 


-S, 


v - k 


dk]ds 


G.  = J 

S. 

J 

The  complex  velocities,  in  (D-l),  are 
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Making  use  of  the  relation  (oC=  slope  of  the  segment) 
dC  = dse ' a 
d£  = dse  'a 

we  evaluate  the  above  integral  with  respect  to  S; 
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Separating  the  above  into  real  and  imaginary  parts, 
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